First results about on-ground calibration of the Silicon Tracker for the AGILE satellite 
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Abstract 

The AGILE scientific instrument has been calibrated with a tagged y-ray beam at the Beam Test Facility (BTF) of the INFN 
Laboratori Nazionali di Frascati (LNF). The goal of the calibration was the measure of the Point Spread Function (PSF) as a 
function of the photon energy and incident angle and the validation of the Monte Carlo (MC) simulation of the silicon tracker 
operation. The calibration setup is described and some preliminary results are presented. 

•Key words: artificial satellites - gamma rays: observations - instrumentation: detectors - telescopes 
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1. The AGILE mission 

AGILE (Astro-rivelatore Gamma a Immagini LEggero) is 
a Small Scientific Mission of the Italian Space Agency (ASI) 
launched on April 2007 and dedicated to high-energy astro- 
physics 111]. The AGILE satellite is designed to detect and im- 
age photons in the 18 - 60 keV, 30 MeV - 50 GeV and 350 keV - 
100 MeV energy bands with excellent spatial resolution, timing 
capability, and large field of view. 

AGILE is the most compact (x 0.25m 3 ), light (120 kg for 
the instrument, 350 kg for the whole satellite) and low power 
(« 60 W) scientific instrument ever developed for high-energy 
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astrophysics. 

The AGILE scientific payload (shown in Fig[T) consists of three 
detectors with independent detection capability. The Gamma- 
Ray Imaging Detector (GRID) consists of a Si-W converter- 
tracker 021] sensitive in the y-ray energy range 30 MeV - 50 
GeV, a shallow (1.5 Xq on-axis) Csl Calorimeter [3] and a seg- 
mented AntiCoincidence system based on plastic scintillators 
fl. 

In addition to the GRID, a coded-mask hard X-ray imaging sys- 
tem (SuperAGILE), made of a Si detector plane and a W mask, 
ensures coverage in the range 1 8 - 60ke V . 
The AGILE main feature is the combination of two co-aligned 
imaging detectors (SuperAGILE and GRID) sensitive in the 
hard X-ray and in the y-ray ranges with large field of view 
l.Osr and « 2. 5sr respectively). 
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Moreover the Csl MiniCalorimeter (MCAL) can operate in 
stand alone "burst mode" in the 350 keV - 100 MeV range to 
detect GRB. 

On ground and subsequently on flight calibrations of a detector 
are essential to the interpretation of its results. The purpose of 
the calibration of a scientific instrument is to reproduce, under 
controlled condition, the detector response in operation. 
This paper describes the on-ground calibration of the silicon 
tracker and some results on the instrument performances de- 
rived by it. 
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Figure 1 : A schematic view of the AGILE scientific instrument. 



2. The Silicon Tracker 

The core of the GRID is the Silicon Tracker (ST) that con- 
verts the y-rays and measures the trajectories of the resulting 
e + /e~ pairs The ST consists of 12 trays with distance be- 

tween middle-planes equal to 1 .9 cm optimized by simulation. 
The first 10 trays consist of a W converter layer 245 fim thick 
followed by pairs of single sided Si microstrip planes with strips 
orthogonal to each other to provide three dimensional points 
(corresponding to a total thickness 0.01(5/) + 0.07(W)X ). The 
last two trays have no W converter layers since the GRID trig- 
ger logic requires at least three contiguous Si planes. 
The detector unit is a 9.5 x 9.5cm 2 tile, 410/im thick with strip 
pitch 121 mum. Four tiles bonded together form a 'ladder'. Ev- 
ery ST plane consists of four ladders. 

Only every second strip is readout to limit the power consump- 
tion. The non readout strips contribute to the resolution through 
the principle of capacitive charge division. 
Each ladder is read-out by three TAA1 ASICs, each operat- 
ing 128 channels at low noise, low power configuration (< 
400^/W/channel), self-triggering ability and analog readout. 
The ST position resolution is below 40^/m for a large range of 
particle incidence angles @]. 



2.1. The GRID simulation 

The GRID as mounted on the spacecraft and as installed in 
the test beam is simulated using the GEANT 3.21 package J3l. 
This package provides for a detailed simulation of the materi- 
als and describes with high precision the passage of particles 
through matter including the production of secondary particles. 
The simulation output is formatted to be readable by the recon- 
struction programs used for the analysis of in-flight data. 

2.2. Direction and Energy Reconstruction 

The y-ray direction reconstruction is obtained from the iden- 
tification and the analysis of the e + /e~ tracks stemming from 
the conversion vertex. Each microstrip silicon plane measures 
separately the X and Y hit coordinates. 

The first step of the event analysis requires to find two tracks 
among the possible associations of the hits detected by the ST 
layers. 

The second step consists in fitting the track trajectories through 
the hits accounting for the presence of energy loss and multiple 
scattering. These steps are performed separately for the X and 
Y coordinates producing four tracks, two for each projection. 
The three dimensional direction is obtained requiring a correct 
association of the two projections of each track. 
The track parameters are fitted by a Kalman filter smooth algo- 
rithm A special implementation of the filter [9] exploits the 
measurement of the angular scattering of the e ± due to the inter- 
actions with the material to estimate the track energies. Com- 
bining the track energies the y-ray energy is estimated. 

3. The y-ray Calibrations 

3.1. Calibration goals 

The goal of the calibration is to estimate the instrument re- 
sponse function by exposing it to a y-ray beam with energy and 
direction known to an accuracy better than the resolving power 
of the instrument. 

The required accuracy of ST is driven by its use during the 
AGILE mission: the systematic errors introduced by the cali- 
bration should be smaller than the statistic errors expected from 
a bright celestial source. 

The detector properties to be evaluated by the calibration are: 
the detection efficiency, the angular resolution, the energy res- 
olution. In this paper we concentrate on evaluating the Point 
Spread Function (PSF) as a function of the y-ray energy and 
incident angle. 

The calibration is also intended to validate the MC simulation 
program. This simulation will be required to complement the 
calibration data in the untested parts of parameter space. In 
particular the information above the maximum energy available 
at BTF can be obtained only through the simulation. 

The calibration is designed to cover a wide range of the ge- 
ometries and conditions realized in space. The ST was cali- 
brated at the INFN LNF in the period 2-20 November 2005, 
thanks to the collaboration between the AGILE Team and 
INFN-LNF. 
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3.2. Calibration strategy 

To meet the calibration accuracy requirements, we have de- 
termined the number of y-rays required for the calibration of 
AGILE, taking into account the photon fluence of a character- 
istic y-rays reference source as Vela. 

With a cover-up efficiency of 50% and an effective area of 
as 500cm 2 , the number of counts estimated is about 10 4 for E 
> 100 MeV, after two months of observation. The requirement 
on the number of calibration photons detected by the GRID is 
4 x 10 5 for E y > 30 MeV, 4 x 10 4 for E y > 100 MeV. 

3.3. Calibration set up 

3.3.1. The Beam Test Facility 

For the ST calibration we used the Beam Test Facility (BTF) 
in the Frascati DAONE collider complex, which includes a 
LINAC at high e~ /e + currents, an accumulator of e~/e + and 
two accumulation rings at 510 MeV. 

The e + je~ beam from the LINAC is led into the accumulation 
ring to be subsequently injected in the principal ring. When the 
beam is not transferred in the accumulator, it can be transported 
from the LINAC in the test beam area through a dedicated trans- 
fer line: the BTF line. The BTF provides a collimated beam of 
e~/e + in the energy range 20-800 MeV with a pulse rate of 50 
Hz. The pulse duration can vary from 1 to 10 ns and the number 
of particles for bunch can range from 1 to 10 5 . 
We operated with energy beam of 463 MeV and a pulse dura- 
tion of 2 ns. 

3.3.2. Target 

y-rays were produced by Bremsstrahlung in a thin Silicon 
target; subsequently a magnet bent away the e~ while the y- 
rays could impinge on the GRID. 

The target is constituted by two pairs of silicon microstrip sin- 
gle sided detectors of 8.75 x 8.75cm 2 and 0.41mm thick, in- 
cluding 384 strips with 228yura pitch. The target measures the 
passage of the <?~ and cause the emission of Bremsstrahlung 
y-ray. 



Figure 2: A schematic of the y-ray line: the target, the bending magnet and the 
PTS. 



3.3.3. Tagging system 

Our Team developed and installed in the BTF area a Photon 
Tagging System (PTS) for the detection of the particles inter- 
acting with the target. The <?~ are tagged using microstrip Si 
detectors located on the internal walls of the bending dipole 
magnet (see Fig(2j). Depending on the energy loss in the tar- 
get, the e~ impinge on different strips. The correlation of the 
measurements of the e~ by the target Si planes and by the PTS 
tags the photon; the position on the PTS measures the photon 
energy. 

The PTS operates in self-trigger mode, i.e. it is readout inde- 
pendently from the GRID. This point has a great relevance for 
the following analysis. 

During the 18 days of calibration, about 2, 10 5 tagged y-rays 
were produced, of which a 40% interacted with the GRID. 

3.4. Instrument Ground Support Equipment ( GSE) 

We developed and installed specific equipment required to 
coordinate and, whenever possible, automate the instrument 
management and the data gathering and analysis as required 
by the calibration procedures. The Mechanical Ground Support 
Equipment (MGSE) 11011 hosts the payload and allows the pre- 
cise motorized translations and manual rotations of the detec- 
tor volume in front of the beam. In near real time, the Science 
Console (SC) 111 111 archives all the instrument data and performs 
the quick look to check the instrument behaviour. It is also in 
charge of producing the energy histogram of the PTS data to 
verify the actual statistic of the PTS measurement and decide 
the measurement duration. 

3.5. Trade-off on the number of e~ /bunch 

The GRID performance should be evaluated in a 'single- 
photon' regime without simultaneous multi-photon interac- 
tions. Multiple photon events are not representative of astro- 
physical conditions and may introduce a significant bias in the 
measurement. 

The best configuration was with 1 e~/bunch, but considering 
the time available for calibration and to obtain a higher effi- 
ciency we adopted 3 e~/bunch. 

3.6. Simulation 

The overall system including the beam terminal section, the 
target, the bending magnet, the PTS and the GRID are simu- 
lated in detail using GEANT 3.21 package [7]. 
That allows a direct comparison between the resolutions mea- 
sured in simulated and real data providing a check of the quality 
of the MC simulations. 

A significant improvements of the comparison between data 
and MC were obtained by overlapping a uniform flux of low 
energy y-rays to the Bremsstrahlung y-ray. These y-ray rep- 
resents a background that cannot be precisely and is tuned to 
match the experimental data. 
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4. Data Analysis 

4.1. Data Samples 

The y-ray beam was directed to the ST at fixed 9 and <p with 
respect to the detector planes. The beam spot on the detector is 
small, « 2 - 3 mm as measured in the target. Yet, ideally the 
photon beam should illuminate uniformly the ST. 
An approximated uniform illumination is obtained by translat- 
ing the GRID with respect to the beam on a run by run basis. 
The beam impinge in four to eight different positions per orien- 
tation, called spills. 

Data were collected for different combinations of 9 = 
0°,30 o ,50° and0 = 0°,45°, 135°, 225°, 270°, 315°. 
Runs for different spills and same orientation are grouped to- 
gether. Also runs for different and same 9 are grouped to- 
gether after having verified that they are compatible. 

4.2. The GRID trigger 

The GRID trigger for AGILE operation is described in [Q]]. 
The relevant point for this calibration is the following: the 
GRID is self-triggering, that is no external signal nor phase 
locking with the accelerator is present. 

During the calibration the in flight trigger configuration was ac- 
tive except for the AntiCoincidence veto that was turned off. 
This choice was imposed by the high rate of background in- 
duced hits in the experimental hall that was reducing the live 
time to an unacceptable level. In this configuration the trig- 
gered events were contaminated by charged particles crossing 
the AntiCoincidence panels. 
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Figure 3: Relation between ST and PTS energy 



4.3. Event reconstruction: the filter 

The event reconstruction program, called filter, selects events 
with a y-ray converting in a e + /e~ pair. The events are expected 
to have two tracks from a vertex within the ST. The kinematic 
of the event is reconstructed applying a Kalman filter. The en- 
ergy of the tracks are estimated by the multiple scattering in the 
ST planes. 

In addition the filter returns a flag assigning an estimation for 
the event being a y-ray or background. 
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Figure 4: E y measured from the GRID fitted with po + p\E P2 for 9 ■ 
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There are four flags tagged as G (gamma), L (limbo), P (par- 
ticle), S (single). The events flagged as G satisfies very strict 
requirements to be a converted y-ray. 
Those flagged as limbo are possible but not certain y-ray. 
Those flagged as particle are estimated to be particle crossing 
the ST (e.g. cosmic muons). 

Those flagged as single are estimated to be single particles from 
a vertex within the ST. 

The Point Spread Function (PSF) of the ST ideally should be 
studied with a sample of G events having a minimal background 
contamination. However, when the PSF is studied for each ori- 
entation and versus the energy, sufficient statistic is required. 
The unselected triggered events at 9 = 0° are more than 3 10 6 
events. The flag fractions are: G(2.3%), L(48.3%) P(38.6%), 
S(10.8%). The striking feature is the low fraction of G events. 
This is a feature of the high background environment present in 
the BTF In this regard the in flight environment is much cleaner 
and the fraction of G events is much higher. 

4.4. Event selection for data: the PTS 

The PTS can be used for two different but related pur- 
poses: 1) as an off-line trigger to identify the emission of 
Bremsstrahlung y-ray in the target and 2) as a device to mea- 
sure the y-ray energy regardless from the ST. 
The PTS and GRID events are paired off-line exploiting the 
event times measured in both devices up to Ips precision. 
The PTS is required to have a very clean signal to reduce multi- 
y-ray events and various background sources. That implies a 
low efficiency as off-line trigger. In the configuration 9 — 0° 
the tagged events are only 23596. The fractions of events in the 
four flags are: G(5.2%), L(44.5%) P(30.3%), S(20.0%). There 
is a significant increase in the fraction of G events that never- 
theless remain a small fraction. The same pattern is present for 
the other orientations. 

The other task of the PTS is the measurement of E y . This is ob- 
tained calibrating with the MC the relation between E 7 and the 
position of interaction of the <?~ on the PTS. A close relation 
between E y measured by the PTS and by the ST is expected. 
Fig[3] shows the profile plot of ST energy versus PTS energy. 
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The correlation is significant, but the spread is large and the lin- 
earity is poor. The GRID energy resolution cannot be evaluated 
precisely with this method. 

The quality of the GRID energy measurement can be estimated 
from Fig|4]where the E 7 spectrum measured by the GRID is fit- 
ted with the function po + as expected for a Bremsstrahlung 
spectrum. The limited distortion indicates that the GRID en- 
ergy resolution and the energy dependence of the efficiency do 
not alter significantly the Bremsstrahlung spectrum. 

4.5. Event selection for data: the phase approach 

A drawback of the PTS approach is the low efficiency « 1%. 
If high quality reconstruction (flag G) is required, the number 
of events available for the PSF determination in any given con- 
figuration may become very small. 

An alternative approach consists in exploiting the BTF bunched 
periodicity at 50 Hz. The intra spill period is Tbtf - 20 ms. 
That implies that subsequent beam related y-ray are spaced in 
time of multiples of Tbtf, in phase with the beam period. The 
event time on the GRID is measured with a resolution of 1 fis 
that defines the precision of the selection. 
In Fig|5]the distribution of time differences between consecu- 
tive events is displayed, showing high peaks in correspondence 
of Tbtf multiples. 

Fig |6] presents the event phase versus the event time. Events 
are in phase if the time difference between consecutive events 
is an integer multiple of Tbtf within lOOyus. These events are 
marked in lighter color. 

Another prominent feature of Fig|6]can be interpreted as fol- 
lows: there are time intervals (approximately 0-500 s and 1800- 
2700 s) showing no accumulation of events in phase. That is 
a sign of beam off time when the GRID measures only beam 
unrelated background. Outside these intervals there is an accu- 
mulation of events in phase with decreasing numbers when the 
number of Tbtf increases. Restricting to the beam on intervals, 
the fraction of events in phase is « 25%. 
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Figure 6: Time difference between consecutive events versus time for a run. 
Events in phase are marked in lighter color 
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Figure 7: PSF (68%) versus energy for MC. 
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Figure 5: Time difference between consecutive events in a run 

Events in phase can be selected regardless the presence of 
tagging to enhance the available statistic and the PSF can be 
estimated by these samples. 

If this approach is correct all tagged events are expected to be in 
phase. That is the case confirming the validity of the approach. 
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Figure 8: PSF (68%) versus energy for data. 
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Figure 9: Gaussian PSF versus energy for MC. 
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the Gaussian PSF should reflect more directly the quality of the 
GRID simulation, rather than the beam simulation. 
An interpretation of these results is that the compatibility of the 
Gaussian PSF for data and MC represents a validation of the 
GRID simulation within the experimental requirements. 

6. Conclusions 

This paper presents some preliminary results of the calibra- 
tion of the AGILE ST at the BTF of the LNF in 2005. 
The setup is described in detail as well as the calibration re- 
quirements. We discussed the problems encountered in exploit- 
ing the PTS originally designed and a novel approach devised 
to circumvent those problems: the phase analysis. 
We concentrated on the measurements of the PSF presenting 
two possible definitions: the Gaussian and the 68% PSF. 
The calibration results are compared with the MC simulations 
for a broad set of variables, showing good consistency with 
some poorer agreement for 68% PSF mainly at low energies. 
These results give confidence on the use of the MC simulation 
in the untested part of the y-ray parameters (e.g. higher E y ) 
especially in flight conditions, i.e. without low energy back- 
ground, and in the measurement of detector parameters, like 
absolute efficiency and energy resolution, that are difficult to 
measure without exploiting the PTS information. 
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Figure 10: Gaussian PSF versus energy for data. 

5. Comparison of PSF for MC and real data in phase 
The PSF was evaluated in two different ways: 

• a Gaussian fit plus a polynomial background, identifying 
the PSF with the Gaussian <x 

• the PSF is identified with the angular spread including 
68% of the events 

The result of the 68% estimation for the PSF for various ver- 
sus E y is shown in Fig|7]for MC and in Figj8]for real data. 
The result of the Gaussian estimation for the PSF for the same 
configurations is shown in Fig|9]for MC and in Fig[l0]for real 
data. 

The 68% PSF is significantly larger than the Gaussian PSF as 
expected in presence of background. 

The data and MC Gaussian PSF are compatible with each oth- 
ers within the statistical errors. 

For the 68% PSF, the data show somehow larger values espe- 
cially at low E y . That is likely due to the low energy y-ray 
background that is not adequately simulated. On the other hand, 
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